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approximants , etc.  A very  considerable  effort  was  expended 
toward  this  end.  The  growing  realization  that  none  of  these 
sophisticated  mathematical  tricks  could  unlock  the  solution 
finally  lead  the  author  to  seek  a new  representation,  ab  initio. 
Furthermore,  it  was  felt  that  the  proper  test  ground  should  be 
the  diatomic  molecule  rather  than  the  triatomic  molecule.  To 
this  end,  the  problem  of  the  vibrational-rotational  interaction 
in  a diatomic  molecule  with  high  angular  momentum  has  been  com- 
pletely re-examined  using  a nonperturbative  scheme  which  completely 
avoid  the  usual  anharmonic  oscillator  interpretation  and  its 
attendant  difficulties.  This  material  is  outlined  in  Section  One. 

Section  Two  is  devoted  to  the  evaluation  of  the  mass- 
dependent  contribution  to  the  potential  energy  in  the  Watson 
molecular  vibration-rotation  Hamiltonian.  An  explicit  solution 
is  obtained  via  generalized  Lucas  polynomials  and  has  been 
checked  against  the  solution  obtained  by  Rothman  and  Clough.  An 
abortive  attempt  to  employ  Chebyschev  operator  expansions  to  the 


van  Vleck  transformation  is  discussed  in  Section  Three. 
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VIBRATION-ROTATION  INTERACTION  IN  DIATOMIC  MOLECULES  WITH 
REFERENCE  TO  LARGE  0 VALUES 
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Schroedinger's  equation  for  the  vibratior.-rotation  inter- 
action of  a diatomic  molecule  is 

-X"  + [Ml+q)-2  + V(q)]X  = AX  (1.1) 

where  A = J(J+1)  and  X(q)  = qi|/(q). 

The  conventional  approach  by  Dunham  [3 ] is  via  perturbation 
theory.  Both  V(q)  and  (l+q)~2  are  expanded  in  a power  series  in 
q.  It  is  possible  to  eliminate  the  constant  and  the  linear  term 
in  q,  so  that  Eq.  1.1  can  be  expressed  as 

-X"  + (aq2  + bq3  + cq4  + • • • )X  ■ AX  (1.2) 

which  is  the  equation  of  the  anharmonic  oscillator.  The  argument 
behind  the  perturbation  approach  is  certainly  plausible,  one 
expects  the  terms  bq3,  cq1*,  •••  to  be  small  compared  to  aq2 . 
Explicit  expressions  for  the  energy  levels  have  been  worked  out 
in  an  elegant  manner  by  Kilpatrick  [2].  However,  these  pertur- 
bation expansions  are  not  convergent! 

In  order  to  see  this  without  undue  complications,  let  us 
confine  ourselves  to  the  case 

-X"  + (q2  + cq4 )X  = AX  . (1.3) 

Straightforward  application  of  the  Rayleigh-Schroedinger 


i 
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perturbation  technique  yields 


A = 1 + if  c - 


c — • • • 


(1.4) 


for  the  ground  state  energy.  The  coefficients  do  not  decrease  in 
magnitude  and,  in  fact,  the  series  does  not  converge  for  any 
nonzero  value  of  c.  Failure  of  the  series  to  converge  follows 
from  the  fact  that  a nonzero  radius  of  convergence  would  have  to 
include  negative  values  of  c ( but , for  negative  c,  the  potential 
energy  does  not  possess  a lower  bound  as  c -*■  -<*>).  The  spectrum, 
consequently,  is  a continuum  and  cannot  connect  continuously 
with  the  discrete  spectrum  for  positive  c.  The  series,  Eq.  1.4, 
is  asymptotically  valid  as  c + 0.  In  view  of  these  remarks,  we 
seek  a non-perturbative  solution. 

An  obvious  approach  is  to  expand  the  wave  function  in  a 
linear  combination  of  harmonic  oscillator  wave  functions;  Chan 
and  Stelman  [3],  Reid  [<#]•  However  (in  part  preparation  for  our 
new  approach),  we  will  assume  that 


X(q)  = £ *«" 

n=0 


(1.5) 


for  the  even  parity  states  and 


00 
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X(q)  = e“q2/2  £ B q2n+1  (1.6) 

n=o  n 

for  the  odd  parity  states. 

Direct  substitution  of  these  series  into  Eq . 1.3  yields  the 
recurrence  relations 

cAn_2  - An_1  + (l+iJn-A)An  - 2n(2n-l)An+1  = 0 (1.7) 

cBn_2  * Bn_x  + (3+1ln-A)Bn  - 2n(2n+l)Bn+1  = 0 (1.8) 

with  A^,  B^  set  equal  to  zero  if  k is  a negative  integer.  Both 
recurrence  relations  are  four-fold  and  the  eigenvalue  problem  is 
the  condition  that  a nontrivial  solution  for  the  (An),  (B^}  exist 
This  leads  to  the  requirement  that  the  A must  be  found  which 

/\  A 

cause  the  vanishing  of  the  infinite  determinants  Mo  and  Mg  of 
which  the  5 x 5 finite  sections  are: 

1 

— 

SEE  SEPARATE  SHEET  (1.9) 


and 


SEE  SEPARATE  SHEET  (1.10) 
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The  determinants  are  banded  and  only  of  width  four  which  is 
an  obvious  factor  in  our  favor.  The  determinant  is  not  symmetric, 
but  this  is  not  a serious  problem  since  there  are  a number  of 
packaged  algorithms  to  evaluate  eigenvalues  of  nonsymmetric 
determinants.  In  order  to  test  the  effectiveness  of  the  approach, 
the  first  four  even  and  the  first  four  odd  eigenvalues  were  cal- 
culated via  12  x 12  finite  sections  of  the  corresponding  infinite 
determinants.  Although  no  serious  attempt  was  made  to  obtain 
very  high  accuracy,  the  results  listed  in  Table  1.1  are  in  rea- 
sonable agreement  with  the  now  standard  results  of  Reid  [4]  based 
on  an  harmonic  oscillator  basis  set  of  one  hundred. 

There  is  no  great  difficulty  in  adding  additional  higher 
order  terms  to  the  potential,  the  main  effect  is  to  make  the  band 
of  the  eigenvalue  determinant  larger.  Consequently  if  one  per- 
sists in  employing  the  anharmonic  oscillator  interpretation,  then 
this  nonperturbative  approach  is  definitely  useful.  The  main 
difficulty,  however,  is  buried  in  the  coefficients  of  the  effec- 
tive potential  since  a,b,c,***  are  functions  of  A.  Consequently, 
large  values  of  A (which  is,  after  all,  the  real  case  of  inter- 
est) mean  that  more  and  more  coefficients  of  the  effective 
potential  must  be  included.  This  in  turn  means  that  the  banded 
matrix  becomes  wider  and  computational  difficulties  in  evaluating 
the  eigenvalues  arise  because  all  the  banded  elements  are  in 
turn  numerically  large. 
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ENERGY  LEVELS  OF  THE  QUARTIC  OSCILLATOR.  I,  LEVELS 
AS  COMPUTED  BY  REID;  II,  LEVELS  AS  COMPUTED  BY 
PRESENT  METHOD. 


TABLE  1.1 


1.060362 

3.799673 

7.455698 


1.060362090^8 
3.79967302980 
7 . 45569793799 
11.6447455114 
16.2618260189 
21.2383729182 
26.5284711837 
32.0985977110 


11.64474 


16.26183 

21.23837 

26.52848 

32.09860 


37.92300 


37.9230010270 


* 


For  this  reason,  we  feel  it  is  expedient  to  abandon  the 
anharmonlc  oscillator  interpretation  and  seek  some  other  repre- 
sentation which  is  more  natural  to  the  intrinsic  problem. 

Let  us  return  to  Eq.  1.1,  rewrite  it  in  the  form 

- ( 1+q ) 2X"  + [A  + ( 1+q) 2V(q) ]X  = A(l+q)2X  (1.11) 

and  deal  with  it  directly.  We  feel  it  is  more  convenient  and 
instructive  to  study  Eq.  1.1  ab  initio  in  the  form  Eq.  1.11  than 
appeal  to  the  anharmonic  oscillator  interpretation.  Our  main 
interest  is  in  large  values  of  A.  To  this  end  we  take  as  our 
prototype  equation 

-(l+q)2X"  + [A  + (1+q ) 2 (bq 2+cq 3+dq4 ) ]X  = A(l+q)2X  . (1.12) 

In  order  to  solve  Eq.  1.12,  assume  that 

00 

X(q)  = e‘q  ' £ C qn  (1.13) 

n=0 

where  the  C are  as  yet  unknown.  Unlike  the  previous  case  of 
n 

the  pure  quartic  oscillator,  the  expansion  must  include  both  odd 
and  even  powers  of  q.  Upon  substituting  the  series  into  Eq. 

1.12,  we  find  (after  some  tedious  manipulations)  that  the 
obey  the  recurrence  relation 
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- f - ■ ~ ' 


g__ 

Y 

■ 


— 


-U+l)U+2)C4+2  - 2t(i+ l)Cl+1  + (-Z2+3l+2+A-A)Cl 


+ (H-2-2A)Cl_1  + (2l-H+a-A)CIL_2  + (0-2)0^ 

+ (Y-DC)l_4  + 6Ct_5  + eC^6  « 0 

(1.14) 

1 

where  C k = 0 for  k = 1,2,***.  Note  that 

ct  = b , y = b+2c+d  , e = d 

0 = 2b+c  3 6 = c+2d  . 

This  is  a nine  term  recurrence  relation. 

The  eigenvalue  conditions  are 

(1.15) 

(M  - AN)C  = 0 

(1.16) 

A 

where  C is  the  column  matrix  of  the  as  yet  unknown  C 

* n 

A A 

. The 

matrices  M and  N are  infinite  matrices.  The  10  x 10 
section  of  M is  given  by 

finite 

SEE  SEPARATE  SHEET 

(1.17) 

while  the  7 x 7 finite  section  of  N is 

1 

SEE  SEPARATE  SHEET 

(1.18) 
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2+A 


(1.17) 


This  is  a somewhat  different  eigenvalue  problem  than  usually 
encountered  in  spectroscopy  in  that  the  eigenvalue  A appears  not 
only  on  the  principal  diagonal  but  also  on  the  first  upper  and 
lower  diagonals.  In  addition,  the  banded  matrix  M is  not  symme- 


tric 


Equation  1.16  is  an  example  of  the  general  eigenvalue  prob- 
lem which  has  been  of  recent  interest  to  numerical  analysis  [5] 
Fortunately  our  problem  can  be  reduced  to  a simpler  eigenvalue 
problem  by  noting  that  the  banded  matrix  fi  is  nonsingular  (i.e.< 
det  N 0),  consequently  its  inverse  N“l  exists.  Premultiplying 
Eq.  1.16  by  N“l  reduces  it  to 


which  is  a regular  eigenvalue  problem  for  a nonsymmetric  matrix 
A = N_1M.  Numerical  calculations  will  be  carried  out  in  the  near 


future 
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2. 


MASS-DEPENDENT  CONTRIBUTION  TO  THE  POTENTIAL  ENERGY  IN  THE 
WATSON  MOLECULAR  VIBRATION-ROTATION  HAMILTONIAN 

Among  the  many  contributions  made  by  Watson  [2]  in  his  im- 


portant paper  on  the  Darling-Dennison-Wilson  Hamiltonian,  Watson 
showed  that  the  effective  reciprocal  inertia  tensor  yag  (a  8 * 1,2,3) 
can  be  expanded  as  a Taylor  series  in  the  normal  coordinates.  The 
use  of  matrix  notation  for  yag  as  a 3 x 3 matrix  y enables  the 
expansion  to  be  cast  into  succinct  form.  The  y matrix  has  the  form 

y = ( 1 0 ) _Js  (E  + | b)”2  (I,)'1*  (2.1) 

where  E is  the  unit  matrix.  Note  that  E and  b commute,  so  we 
can  use  the  binomial  expansion  of  (E  + 1/2  b)-2;  thus 


y = d#)  14 

E + 2 Crtl)  (-i)r  br 

(I#)_ls  . (2.2) 

r-1 

A 

The  matrix  b is  also  3 x 3,  consequently  we  can  express  any  power 
of  b greater  than  or  equal  to  three  in  terms  of  the  generalized 
Lucas  polynomials  of  order  three  introduced  by  Barakat  and 
Bauman  [2]  to  handle  just  such  problems. 

The  generalized  Lucas  polynomials  of  order  three  are  defined 
by  the  recurrence  relation  [2] 


( 3 

) _ $ u\,(*  , d>  , 

y2  n+1  y 


♦.) 


(2.3) 
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subject  to  the  initial  conditions 

d!')(*i»*i»*»)  * 0 » - 1 


D»<,)(VW  = 0 , UJ(,)(*1, ^ . 


(2.4) 


The  explicit  solution  of  Eq.  2.3  subject  to  these  initial  condi- 
tions is 


1 » J 

where  i + 2j  + 3k  = n.  The  <j>'s  are  the  symmetric  functions 
formed  by  the  traces  of  6,  b2,  b3.  The  traces  are  invariant 
under  a unitary  transform  and  are  thus  independent  of  orientation. 
The  ♦'s  are  related  to  the  o’s  (the  traces  of  powers  of  b,  i.e., 
ok  = trfck)  by 

♦i  “ °i 


2 ^ai 


f 


3oioj  + 2os) 


(2.6) 


Now,  if  we  set 


A A 

b -2wb.  Q. 


(2.7) 


where  {Qk>  are  the  normal  coordinates,  then 
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bXX 

•<? 

. xz 
bk 

s,.- 

byy 

byz  . 

k 

(2.8) 

hzx 

bk 

KZZ 

bk 

Straightforward  calculations 

then 

yield 

,,h’1  -?ek‘ 


tr  b2  = o q 

tr  D " °2  k l gkJt  wkw* 


tr  £3  = 0j  =££Eg(^QkQAQm 


where 


(2.9) 

(2.10) 

(2.11) 


kxx  x >,yy  x kzz 
bk  + bk  + bk 


.xx.  xx  Kyy»,yy  . Kzz.  zz  x o^y^y2 

bk  b*  + bk  V + bk  b*  + 2bk  b* 


ovxz.  zx  ,hyzKzy 
+ 2bk  b^  + 2bk  b^ 


xx  xx  xx  byybyybyy  + bzzbzzbzz 
k l m k l m k l m 


+ + 3bjxbfb^x  ♦ 3byybyxbxy 


♦ 3byybyzbzy  ♦ 3bzzbfxbxz  ♦ 3b£zbzybyz 


(2.12) 


(2.13) 


k ul  m 


k “1  m 


x ovXywyzKzx  . ,KxzKzy.  j 
+ 3bk  bl  bm  + 3bk  bt  br 


(2.1M) 


Consequently,  ♦ j,  $2,  and  <frs  can  be  expressed  directly  In  terms 
of  the  {Qk>,  the  explicit  formulae  are 


o,  = tr  b 


,(l)Q 

k wk 


(2.15) 


o2  = tr  b2 


E hkk)QkQi 
k,£  KK  K *• 


(2.16) 


’>  = tr  ba  - £ 

k,£,m 


(2.17) 


with 


.(»)  = _(») 


h(2>  = i . i g<2> 

nkl  2 °k  g£  2 gk£ 

h(.3)  l gOy.y.j  . 1 .(»)  (.)  + i .(>) 

ktm  Z 6k  g£  ®m  2 gk  g£m  3 sk£m 


(2.18) 


(2.19) 


(2.20) 


Note  that  ai  is  linearly  dependent  on  {Q^},  o2  is  a quadratic 
function  of  (Qk>,  while  a3  is  a cubic  function  of  (Qk).  All 
higher  o's  can  be  written  in  terms  of  , o2,  and  o3.  In  fact. 


we  can  prove 


°k  ■ “i”0. + * 3*.u£ii 


(2.21) 


for  k > 3. 
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Barakat  and  Bauman  [2]  have  shown  that  any  positive  power  of 
6 greater  than  or  equal  to  three  can  be  expressed  as 


br  = u£3)b2  + [u£3j  - ♦1U^*)]b  + *SU^E  . (2.22) 


This  equation  is  the  key  to  the  subsequent  analysis.  Let  us  re- 
turn to  Eq.  2.2  and  employ  Eq.  2.22  in  the  right  hand  side,  the 
result  is 


A A 


U = (I0) 


-h 


|E-b  + | b2  + X)  (r+1)  u£3)b3 


00 

♦ E 

r=3 

(r+1) 

(4f 

u(3)  - <f>  u(3 

r+1  M r 

00  „ \ 

♦ E 

r=3 

(r+1) 

(-D 

-h 


(2.23) 


or 


; = (i,)-  - Ci,rhbCi,)-h  * £ a,)-\‘(i,rh  ♦ <r*i>  (-if 

■ ♦ t3  (r*l>  (-if  [u<;>  - ♦,„<■>) 

• (i,)-* bd,rh  ♦ E3  <r+l)  (-|f  <i,)“  ■ <2-2«) 


Actually  we  really  want  to  work  with  the  S-matrices  given  in 


Watson  [1] 


I 


A A _LA  /V 

b = (I)  Vi,)  * 
0 0 


(2.25) 


A _ A A — X>.  A .A  A . 

(I.)  ’b ( I ) * « (I  ) a(I  )-1  . 
0 0 0 0 


(2.26) 


We  can  also  show  that 


(I.)  ^*(1.)  * = d0)-1a(I0)-1a(Io)-1  . 


(2.27) 


Upon  employing  the  expressions  in  Eq.  2.24,  we  can  express  p 
entirely  in  terms  of  the  a.  The  final  result  is 


v = Ad,)-1  + B(I0)-1a(I0)-1  + C(I0)-1a(I0)-1a(I0)-1  (2.28) 


where 


00  r 

A = 1 + 2 B+l)  (*|)  Mr-l  (♦1.*a.*3) 


(2.29) 


00  r 

B = -1  + S (r+1)  (-|)  C*i. ♦*.♦,)  - 


( <l>  j > <J>  2 » <l>  3 ) 1 


0 s i + S, <r+1)  Bf « 


^3)  (♦1»*2»*j)  ‘ 


(2.30) 


(2.31) 


Although  Eq.  2.28  is  the  desired  result,  we  still  have  to 
express  A,  B,  and  C in  terms  of  {Q^}.  This  is  easily  accomplished, 
albeit  tedious.  In  C,  we  encounter  u£3^  alone.  Upon  employing 
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the  explicit  expressions  for  listed  in  Barakat  and  Bauman 

we  can  prove  that 

u£3)  -E  ft2!..k  Qk  Qk  •••  Qk  (2. 

k,  k,  k Kr  K1  Kr 

12  r 


with 


The  method  is  ideally  suited  to  automatic  evaluation  on  a 


computer  and  there  Is  no  need  to  carry  out  higher  order  explicit 
results  by  hand. 

. The  series  in  A contains  the  function 

Wr+2  5 Q3U£-1  » r = 3, (2.39) 
while  the  series  in  B contains 

Vr+1  E Uiil  - *iur3)  * r * 3,4, •••  (2.40) 

The  first  few  members  of  both  sets  have  been  worked  out  explicitly 
and  are 


(2.41) 


(2.42) 


(2.43) 


(2.44) 


(2.45) 
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To  summarize,  the  explicit  forms  of  A,  B,  C are: 


(2.M6) 


(2.47) 


(2.48) 


(2.49) 


(2.50) 

B * I1  ♦ I - tk  V<!')  + ts  vl!>  - & v‘,)  ♦ ik  vl°]  <2-51) 
c - «<*> . £ »!•’  - i *<"  ♦ A - & <3)  + 5?s 


(2.52) 


to  sixth-order  in  the  Q. 


These  results  are  essentially  equivalent  to  those  obtained 
by  Rothman  [3]  using  an  iterative  scheme.  Dr.  Rothman  kindly 
lent  me  some  numerical  results  obtained  via  his  approach,  the 
present  more  explicit  method  yielded  the  same  numerical  results. 


NOTE  ON  OTHER  SCHEMES 


The  author  had  originally  hoped  that  the  method  of  matrix 
PadS  approximants  [4]  could  be  utilized  directly  to  handle  Eq. 
2.1.  Unfortunately  this  approach  is  almost  useless  because  one 
cannot  explicitly  calculate  anything  except  the  first  matrix 
Pad£  approximant.  In  view  of  this  fact,  we  omit  any  discussion 
of  the  method  and  of  our  fruitless  manipulations  directed  toward 
this  end. 
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3.  NOTE  ON  A CHEBYSCHEV  OPERATOR  EXPANSION  OF  THE  van  VLECK 

TRANSFORMATION 

The  van  Vleck  transformation  has  been  applied  to  the  problem 
of  vibration-rotation  interaction  of  polyatomic  (and  particularly 
triatomic)  molecules  by  various  investigators,  for  one  of  the 
more  successful  attempts  see  Rothman  and  Clough  [I]. 

I have  made  an  attempt  to  seek  a different  type  of  expansion 
of  the  van  Vleck  transformation  that  would  converge  more  rapidly 
than  the  usual  power  series  procedure  and  yet  hopefully  retain 
the  versatile  attributes  of  the  power  series.  A method  that 
appeared  to  have  both  properties  is  a Chebyschev  expansion  ex- 
tended to  operators.  It  is  well  known  among  numerical  analysts 
that  Chebyschev  expansions  are  more  rapidly  convergent  than  the 
corresponding  power  series  [2]. 


The  van  Vleck  transformation  is 


H' 


eiXS  g e-iXS 


(3-D 


where  the  operator  S is  hermitian  and  X is  an  expansion  parameter 
whose  value  lies  between  zero  and  unity.  For  our  purposes,  it 
is  more  convenient  to  deal  with  the  "real"  and  "imaginary"  parts 
of  the  right  hand  side  of  Eq.  3.1;  thus 

R = (cosXS)H(cosXS)  + (sinXS)H(sinXS)  (3-2) 

I = (sinXS)H(cosXS)  - (cosXS)H(sinXS)  . (3.3) 
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sinX§  - 2 £ (-1)nj2n+l(X)T2n+l(§) 

n*l 

a 

where  Tn(S)  is  the  Chebyschev  polynomial  of  degree  n in  the  oper- 
ator §.  In  order  that  this  method  work  it  is  manditory  that  the 
norm  of  § be  less  than  or  equal  to  unity,  which  is  a severe 
restriction. 

Although  these  expressions  converge  rapidly  they  cannot  be 
utilized  in  their  present  form  since  they  must  be  rearranged  in 
power  of  X in  order  to  be  useful.  This  requires  that  we  write 
out  the  power  series  expansions  of  Jn(X)  and  collect  equal 
powers  of  X. 

Unfortunately,  when  this  extremely  tedious  procedure  was 
carried  out,  the  resultant  expressions  are  so  involved  that  any 
attempt  to  actually  perform  the  van  Vleck  transformation  was 
frustrated  by  the  cumbersome  formulae.  This  possible  approach 
does  not  offer  any  advantages  over  the  usual  power  series  approach 
to  the  van  Vleck  transformation. 
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